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智能 排 班 问题 的 概率 约束 规划 模型 与 有 效 求解 


陈 志 平 :， 刘 w, OP 
(1- 西安 交通 大 学 理学 院 科学 计算 与 应 用 软件 系 ， 西 安 710049; 
2- 西安 交通 大 学 金条 经济 研究 中 心 ， 西 安 710049) 

摘 要 : 由 于 其 广泛 的 应 用 领域 与 难 解 性 ， 对 不 同形 式 智 能 排 班 问题 的 研究 现 已 成 为 运筹 学 界 的 一 个 研究 
热点 。 为 克服 现 有 模型 无 法 恰当 反映 人 力 需 求 约束 实际 特征 与 求解 算法 效率 低 等 不 足 ， 本 文 对 综 
合 排 班 问题 构建 了 两 种 不 同形 式 的 概率 约束 随机 规划 模型 。 利 用 样本 平均 近似 方法 和 “大 M 方 
法 ”， 我 们 将 概率 约束 问题 转化 为 大 规模 整数 线性 规划 问题 ， 并 综合 运用 了 Lagrange 松弛 法 、 
分 解 算法 与 信赖 域 算法 等 方法 的 思想 ， 设 计 出 了 一 有 效 的 嵌 套 分 解 算法 。 最 后 ， 通 过 数值 算 例 说 
明了 本 文 模型 的 合理 性 与 所 设计 算法 的 有 效 性 。 
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1 引言 


所 谓 排 班 问题 ， 即 是 在 给 定 的 一 系列 要 求 ( 比 如 每 天 的 用 工人 员 需 求 、 工 作 强 度 要 求 、 劳 动 
保障 要 求 等 ) 之 下 ， 对 给 定数 量 的 工作 班组 进行 合理 的 工 休 安 排 。 在 过 去 ， 排 班 往往 是 一 个 人 
工 统 筹 的 过 程 ， 但 随 着 大 型 呼叫 中 心 等 现代 服务 模式 的 出 现 ， 如 此 大 规模 的 排 班 是 手工 不 可 能 
完成 的 任务 。 这 时 ， 利 用 计算 机 找到 一 种 新 的 智能 化 方法 就 成 为 该 问题 的 解决 之 道 。 智 能 排 班 
旨 在 利用 计算 机 与 现代 优化 技术 解决 劳动 力 管 理 问题 ， 为 客户 安排 最 佳 任务 ， 使 其 在 满足 未 来 
需求 预测 的 情况 下 ， 还 能 满足 运营 商 的 各 种 人 性 化 需求 (包括 班组 间 的 均衡 )。 它 在 现代 服务 行 
业 ， 诸 如 发 电 、 电 信 、 人 金融 等 领域 正在 被 广泛 采用 。 对 于 劳动 密集 型 行业 ， 其 运营 成 本 和 运营 
效率 很 大 程度 上 取决 于 如 何 更 有 效 地 安排 员工 的 工作 任务 ， 智 能 排 班 方法 正 是 根据 预测 估计 相 
应 的 人 力 需 求 ， 合 理 安排 任务 ， 以 达到 人 力 资源 最 有 效 利用 的 目的 。 

自 文献 四 提出 道路 收费 亭 问题 的 集合 覆盖 模型 以 来 ， 数 学 规划 技术 已 被 广泛 应 用 于 排 班 问 
题 的 建 模 与 求解 中 。 基 于 对 员工 排 班 的 不 同 刻 画 ， 现 有 的 模型 可 分 为 三 类 : 工 休 排 班 、 班 次 排 
班 和 综合 排 班 加。 工 休 排 班 是 在 一 个 周期 内 为 员工 安排 工作 日 和 休息 日 如， 班次 排 班 仅 考虑 一 
天 的 情况 ， 在 该 天 内 为 员工 安排 上 班 时 间 和 休息 时 间 和 内， 综合 排 班 则 将 工 休 排 班 和 班次 排 班 相 
结合 ， 既 确定 一 个 周期 内 的 工作 日 和 休息 日 ， 又 为 每 个 工作 日 安排 上 班 时 间 和 休息 时 间 句 。 

由 于 其 简单 、 直 观 、 易 操作 等 特性 ， 确 定性 数学 规划 模型 在 员工 排 班 中 占据 研究 与 应 用 的 
主体 。 几 乎 所 有 这 些 模型 均 是 在 最 佳 满足 业务 需求 的 条 件 下 ， 最 小 化 人 力 成 本 。 根 据 具体 的 应 
用 环境 ， 此 模型 有 多 种 不 同 角度 的 推广 ， 典 型 的 包括 更 精细 地 划分 时 间 段 全 、 更 人 性 化 地 分 配 
员工 工作 时 间 团 、 抑 或 加 入 各 种 连续 上 班 和 连 班 间隔 等 约束 回 。 

至 于 如 何 具体 表达 各 种 业务 需求 ， 且 前 主要 有 如 下 两 种 描述 方式 。 一 种 是 利用 单 边 不 等 式 
或 双边 不 等 式 ， 保 证 各 时 段 所 安排 人 数 大 于 某 预 定 值 (如 入 力 需求 的 预测 值 )， 或 限定 其 处 于 某 
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个 范围 内 。 单 边 不 等 式 法 能 有 效 地 保证 任意 时 刻 的 人 力 需 求 都 能 得 以 满足 ， 但 缺点 是 不 能 限定 
每 个 时 刻 应 安排 人 力 的 上 限 ， 这 有 时 会 导致 人 员 宛 余 严 重 ， 而 双边 不 等 式 法 在 单 边 不 等 式 的 基 
础 上 ， 对 每 时 刻 的 人 力 增 加 一 个 上 界 ， 从 而 限制 了 人 员 宛 余 的 程度 ， 但 如 何 选取 这 样 的 上 界 是 
较为 困难 的 问题 。 如 果 选 取得 太 小 ， 则 可 能 导致 模型 无 解 ， 若 选取 得 过 大 ， 又 难以 达到 将 人 员 
宛 余 平 摊 的 目的 。 另 一 种 是 文献 [2,9] 等 中 提出 的 目标 规划 方法 。 这 种 方法 最 大 的 特点 在 于 放松 
了 严格 的 上 、 下 界 要 求 ， 将 实际 安排 人 力 与 预测 人 力 的 偏离 程度 作为 目标 来 最 小 化 。 

虽然 确定 性 模型 有 多 种 人 力 需 求 约束 的 刻画 方式 ， 但 其 确定 性 需求 量 的 假设 严重 限制 了 
它 的 应 用 范围 与 应 用 效果 。 熟 知 ， 智 能 排 班 中 最 核心 的 内 容 就 是 要 最 优 地 满足 不 同 的 资源 需 
求 。 然 而 ， 未 来 的 资源 需求 量 是 事先 无 法 预知 的 随机 量 ， 而 预测 其 最 食 适 的 取 值 又 是 很 困难 
的 89。 鉴 于 此 ， 一 个 自然 的 想法 就 是 直接 将 资源 需求 作为 随机 变量 ， 用 随机 规划 方法 来 描述 
需求 约束 ， 己 有 文献 中 有 两 种 常见 的 处 理 方法 : 一 是 将 人 力 需 求 约束 不 满足 程度 的 期 望 值 作 
为 惩罚 加 入 目标 函数 ， 此 即 随机 目标 规划 法 0; 二 是 将 人 力 需求 约束 满足 的 概率 作为 关注 重 
点 ， 对 满足 的 概率 值 加 以 约束 ， 即 概率 约束 规划 法 213 。 这 两 种 方法 各 有 优 缺 点 。 随 机 目标 
规划 从 平均 的 角度 去 度量 人 力 需 求 约束 的 满足 情况 ， 不 能 有 效 刻 画 随机 需求 分 布 的 尾部 特性 对 
约束 的 影响 ， 概 率 约束 规划 将 不 好 的 事件 发 生 的 概率 限定 在 某 个 范围 之 内 ， 这 样 能 够 根据 实际 
要 求 灵 活 选择 概率 疹 值 ， 但 如 何 设置 合适 的 置信 度 比 较 困 难 ， 且 这 类 约束 没有 考虑 在 约束 限定 
之 外 的 极端 情况 。 

对 所 建立 随机 模型 的 求解 ， 与 现 有 求解 概率 约束 问题 的 方法 一 样 ， 通 常 的 做 法 是 从 随机 变 
量 的 分 布 特性 入 手 ， 将 概率 约束 转化 为 确定 性 约束 ， 再 利用 求解 确定 性 规划 问题 的 相关 算法 予 
以 处 理 。 至 于 如 何 将 概率 约束 确定 化 ， 主 要 有 三 种 思路 : 一 是 在 假设 正 态 分 布 、t 分 布 等 特殊 
分 布 之 下 ， 直 接 将 概率 约束 转化 为 等 价 的 确定 性 约束 84， 二 是 利用 平均 话 务 率 等 参数 ， 基 于 
排队 论 中 采用 的 分 布 模型 来 推导 出 人 力 需 求 量 服 从 的 分 布 b142 ， 再 通过 抽样 方法 将 概率 约束 
确定 化 ， 三 是 直接 假设 随机 人 力 需 求 服从 某 个 离散 分 布 h5161， 将 一 个 周期 内 的 人 力 需 求 模拟 
为 数 种 不 同 的 需求 向 量 ， 每 组 向 量 代表 需求 的 一 组 实现 ， 再 对 每 个 向 量 赋予 一 定 的 概率 (如 等 
概率 )， 使 总 概率 为 1。 以 文献 上 7] 中 的 构造 法 为 例 : 它 将 需求 分 为 周期 内 每 日 需求 模式 与 每 天 
内 每 小 时 需求 模式 两 层 ， 每 一 层 又 可 根据 增长 趋势 、 四 趋势 等 不 同类 型 构造 6 种 模式 ， 再 将 这 
两 层 结合 起 来 产生 离散 分 布 的 样本 。 有 了 离散 分 布 ， 就 可 直接 将 概率 约束 转化 为 通常 的 确定 性 
约束 。 

从 实际 调研 中 我 们 发 现 ， 对 于 现实 中 的 排 班 问题 ， 相 对 于 人 力 需 求 约束 满足 程度 偏差 之 大 
小 ， 客 户 更 为 关注 的 是 偏差 存在 的 情形 是 否 过 多 、 能 和 否 在 大 多 数 情况 下 使 人 力 需 求 得 到 满足 和 
如 何 恰当 描述 并 与 灵活 控制 人 力 需求 约束 等 随机 性 约束 的 满足 。 另 外 ， 对 于 人 力 需 求 等 约束 ， 
硬性 的 下 界 约 束 很 容易 导致 人 员 宛 余 严 重 ， 如 何在 某 些 时 段 的 人 力 不 足 与 另 一 些 时 段 的 人 员 元 
余 间 建立 平衡 以 使 人 力 成 本 最 低 ， 至 今 还 没有 好 的 解决 方法 。 最 后 ， 对 于 现实 中 复杂 的 大 规模 
排 班 问题 ， 现 有 求解 随机 模型 的 算法 效率 比较 低 ， 怎 样 设计 可 快速 求解 的 有 效 算法 是 隔 待 解决 
的 难题 。 

针对 已 有 文献 的 上 述 不 足 ， 本 文 将 就 综合 排 班 问题 的 建 模 与 求解 进行 研究 。 为 灵活 而 恰当 
地 描述 人 力 需 求 约束 ， 我 们 引入 人 力 不 足 量 的 概念 ， 通 过 对 其 设 定 下 界 ， 并 选取 适当 的 置信 水 
平 ， 建 立 两 类 不 同形 式 的 概率 约束 条 件 来 刻画 人 力 需求 约束 ;为 权衡 人 力 不 足 与 人 力 元 余 等 问 
题 ， 我 们 在 建 模 时 引进 加 权 惩 罚 的 思想 ， 不 同 于 现 有 方法 ， 为 有 效 求解 所 建立 的 概率 约束 规划 
问题 ， 我 们 采用 随机 规划 领域 中 新 近 提 出 的 样本 平均 近似 方法 来 将 概率 约束 规划 问题 转化 为 确 
定性 整数 规划 问题 ， 并 利用 “大 M 方法 ”的 思想 将 其 线性 化 ， 为 快速 求解 所 导出 的 大 规模 整数 
规划 问题 ， 我 们 综合 运用 了 Lagrange 松弛 法 、 分 解 算法 与 信赖 域 算法 等 方法 的 思想 ， 设 计 出 
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了 一 有 效 的 嵌 套 分 解 算 法 。 最 后 通过 数值 算 例 说 明 本 文 模型 的 合理 性 与 所 设计 算法 的 有 效 性 。 


2 ”模型 介绍 


为 灵活 反映 实际 中 对 人 力 和 需求 约束 的 不 同 要 求 ， 我 们 将 依据 该 类 约束 所 限定 的 范围 ， 采 用 
两 种 概率 约束 形式 进行 建 模 : 一 种 是 分 别 对 每 个 时 段 加 以 限制 ， 另 一 种 则 是 对 整个 周期 引入 统 
一 的 一 个 约束 。 

考虑 有 了 个 时 段 ，J7 个 班次 的 排 班 问题 ， 我 们 首先 探讨 对 各 个 时 段 分 别 加 以 限制 的 模 
型 。 为 精细 刻画 不 同时 段 上 人 力 需 求 的 满足 程度 ， 将 每 个 时 段 上 所 安排 人 数 与 该 时 段 实 际 
ARKAA r t= 1,… ,全 的 差 定义 为 人 力 不 足 量 ， 同 时 给 每 个 人 力 不 足 量 一 个 下 限 -6, t = 
1,… ,全 。 在 每 个 时 段 上 上 ， 根 据 问题 的 具体 要 求 ， 保 证 人 力 不 足 量 不 低 于 -6 的 概率 不 小 于 某 
个 置信 度 ae 。 以 极 小 化 总 的 人 力 需 求 为 目标 ， 我 们 可 建立 下 列 概率 约束 规划 模型 


J 
min Dz; (1) 
j=l 
J 
s.t. P(Y ajax; -re > 6) > ar, t=1,2,---,T, (2) 
j=l 
rex, (3) 


其 中 7 = 1,… ,J 了 表示 不 同 的 班次 ， 要 求 任意 两 个 班次 均 不 相同 。at; € {0,1} 表示 班次 了 在 时 
段 + 是 否 为 上 班 状态 ，1 表示 上 班 ，0 表 示 不 上 班 。z7 为 决策 变量 ， 即 上 班次 7 MATH. XR 
示 除 人 力 需 求 约束 (2) 外 的 其 它 所 有 约束 。 例 如 ，X 可 包含 如 下 三 种 典型 的 限制 。 

一 是 变量 取 整 及 其 大 于 等 于 零 的 约束 


aj>0, 2j€Z, j=1,2, J. 


二 是 公平 性 约束 ， 即 要 求 周期 内 所 有 上 班 员工 的 平均 工作 时 长 为 预先 给 定 的 值 A 五 
J J 
(S045) = AH x}, 
j=l j=1 


其 中 六 表示 班次 ?在 排 班 周 期 内 的 上 班 时 长 。 
三 是 特殊 班次 限制 ， 即 要 求 不 同班 次 集合 r 中 安排 的 员工 人 数 不 小 于 某 预 先 给 定 的 下 
界 SRPim， 同 时 不 大 于 某 预 先 给 定 的 上 界 S Rx 


Sap ORE Dy cy SSR VIER, 
JEG, JEG, 


其 中 G, 表示 班次 集合 ~ 中 包含 的 班次 。A 表示 所 有 需 满足 相应 类 型 约束 的 班次 集合 的 集合 。 
当然 ， 我 们 还 可 根据 实际 需要 ， 增 加 更 多 更 详细 的 约束 ， 因 这 部 分 约束 并 不 牵扯 随机 因 
素 ， 改 而 可 参考 文献 [18,19] 等 中 确定 性 排 班 问题 的 建 模 方法 。 因 非 本 文 考虑 的 重点 ， 以 下 不 再 
对 其 进行 详细 论述 。 
若 不 对 不 同时 段 赋予 特定 的 满足 概率 值 ， 则 可 采用 统一 的 置信 水 平 ， 并 给 出 另 一 种 形式 的 
概率 约束 表述 : 设 要 求 各 个 时 段 上 人 力 不 足 量 均 不 低 于 最 低 限度 的 概率 不 小 于 某 个 置信 度 a, 
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则 可 构建 如 下 的 联合 概率 约束 模型 
J 
min 》 3; (4) 
jal 
J 
s.t. P( > ayz; — re +ð > 0, t=1,---,T) 2a, (5) 
j=l 
rEX. (6) 


上 述 两 种 模型 (1)-(3) 和 (4)-(6) 与 已 有 文献 中 的 模型 233 相 比 更 加 灵活 ， 可 根据 对 不 同 
天 、 不 局 时 段 的 关注 程度 加 以 不 同 的 偏差 容忍 程度 和 概率 约束 满足 程度 。 同 时 ， 在 入 中 也 综合 
考虑 了 确定 性 排 班 文献 中 较 新 的 成 果 hsa9)。 


3 ”模型 求解 


相 较 概率 约束 规划 模型 (1)-(3)， 联 合 概率 约束 (5) 中 的 概率 事件 并 非 由 单一 线性 不 等 式 所 表 
示 ， 而 是 由 人 个 线性 不 等 式 联 合 表示 。 对 些 ， 可 通过 无 穷 范 数 将 其 表示 为 由 一 个 不 等 式 表 示 的 
概率 事件 的 形式 。 为 书写 方便 并 用 统一 方式 表达 这 两 种 概率 约束 模型 ， 我 们 引入 以 下 记号 。 首 


先 ， 念 
J 
T= (£1, E2 zyj) ER, r= (rT, rr) € RT, a 
j=l 
进而 ， 记 
J 
G(x, Ti) = 》 aux; —ret+, t=1,2,---,T, 
j=l 
q(x) = P(Gi(2, ri) > 0) =E(doo)Gr(z,r1)), t= 1,2,---,T. 
这 里 To,coGe(z,m) 表示 指示 函数 ， 即 当 Gi(z,mr) > 0 时 ，1o Gai(zri) = 1, “4Gi(a,rt) < 
0 时 ，1o coGt(zre) = 0。 类 似 地 ， 记 


J 
G(z,r) = min{ X aijt i + Ôt, t = 1, ets, & Th 
j=1 


q(x) = P(G(a,r) > 0), 
则 概率 约束 模型 (1)-(3) 就 可 转化 为 


min f(x) (7) 
st. “gay > a t=1,2,..,T, (8) 
rEX, (9) 
而 概率 约束 模型 (4)-(6) 则 可 转化 为 
min f(z) (10) 
s.t. g(x) >a, (11) 


TEX. (12) 
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因 转 化 后 的 两 个 模型 具有 相同 的 形式 ， 以 下 仅 就 概率 约束 模型 (7)-(9) 的 求解 展开 论述 。 

为 克服 已 有 算法 对 分 布 类 型 的 依赖 性 、 难 于 实现 的 不 足 ， 我 们 采用 随机 规划 领域 中 新 近 
提出 的 样本 平均 近似 i1429 策略 来 转化 概率 约束 ， 其 最 大 的 优点 在 于 不 需要 假定 原始 的 分 布 类 
型 ， 仅 仅 通 过 抽取 一 组 样本 来 近似 原 随 机 分 布 ， 且 可 保证 收敛 性 。 

对 随机 变量 re， 设 ri,i = 1,… ,TJ， 为 ri 的 I 个 独立 同 分 布 的 样本 。 我 们 用 这 了 个 样本 来 近 
似 其 本 身 的 分 布 ， 并 用 这 些 样 本 的 平均 来 近似 其 原本 分 布 的 期 望 ， 即 用 


I 

a 1 i 

G(x) = 》 710) (Ge(Z, ri) 
i=l 


来 近似 代替 q(z)。 那 么 ， 对 于 给 定 的 概率 或 置信 水 平 Y e (0,1), t= 1,---, 7, E (7)-(9) 则 
可 转化 为 下 列 最 优化 问题 


min f(x) (13) 
st. @(z)>%, t=1,2,---,T, (14) 
ce Xx. (15) 


我 们 将 此 问题 称 之 为 原 问题 (7)-(9) EEKE w, t= 1,… ;三 下 的 样本 平均 近似 问题 。 关 于 
上 述 样本 平均 近似 方法 的 收敛 性 ， 由 文献 [21] 中 的 分 析 ， 我 们 有 如 下 结论 。 

定理 1 假设 原 问题 与 样本 平均 近似 问题 的 置信 度 相同 即 ， =, t=1,---,.T, AR 
概率 约束 外 的 可 行 域 为 紧 集 ，G(z,r) A Caratheodory 函数 ， 且 有 如 王 条 件 成 立 : 设 原 问 题 
KRENT, HEAKKE > 0， 存 在 ze X, EB |r- 2 之 ce 和 gi(x) >a, Mi 一 wv*， 
且 当 了 一 +ooltt, MERIA D(S1,S) > 0, K BG, Sr 分 别 表示 样本 平均 近似 问题 (13)- 
(15) 的 最 优 值 和 可 行 集 ， 而 ur*， 5 分别 表示 原 问 题 (1)-(4) 的 最 优 值 和 可 行 集 ，D(S1, 5) 表示 集 
A Si 与 5 之 间 的 距离 。 

对 比 样本 平均 近似 问题 与 通常 离散 分 布下 所 得 问题 可 以 看 出 ， 样 本 平均 近似 也 可 被 看 做 一 
种 有 具有 了 个 相等 概率 样本 的 特殊 分 布 。 在 置信 度 不 变 的 情况 下 ， 当 了 趋 于 无 穷 大 的 时 候 ， 样 本 
平均 近似 问题 的 解 收敛 于 原 问 题 的 解 、 在 了 足够 大 的 情况 下 ， 若 关 > a, t=1,---,T, WE 
本 平均 近似 问题 的 可 行 解 亦 是 原 问题 的 可 行 解 ， 而 当 < ai; t = 1,… ,了 时 ， 样 本 平均 近似 
问题 的 最 优 值 是 原 问 题 的 下 界 。 所 以 样本 平均 近似 方法 也 可 用 于 构造 可 行 解 或 是 估计 最 优 解 的 
下 界 。 

Al @ (2) 中 含有 非 光滑 的 示 性 函数 ， 导 致 问题 (13)-(15) 难于 处 理 。 为 此 ， 我 们 使 用 类 似 文 
献 [20] 中 的 “大 M 方 法 ”， 通 过 引入 辅助 变量 来 将 示 性 函数 转化 为 线性 形式 。 具 体 地 ， 令 


(16) 
WU Ge(w) 就 可 以 用 zi 的 线性 函数 表示 如 下 
I 1. 
d(x) = S054 
wl 


RE i 5 ox, ri Ze) E REPER Gi (a, ri) + M(1 i) > 0 与 Gi(z,7i) < Mzi, zi € 
{0,1}, HP M 是 一 个 取 值 比 Gi(zx, rt) 的 最 大 值 还 要 大 的 数 。 因 为 在 概率 约束 中 我 们 已 对 愉 加 
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入 了 更 严格 的 下 界限 制 ， 不 等 式 Gi(x,7i) < M 闪 实际 上 已 不 起 作用 ， 所 以 在 模型 中 可 不 加 入 
此 不 等 式 。 

通过 (16) 式 中 的 引入 和 上 述 的 转换 ， 我 们 就 可 将 原 概 率 约束 模型 (1)-(4) 转化 为 如 下 形式 
的 整数 线性 规划 问题 


min f(z) (17) 


I 
1. 
s.t. 2 $= 1,2, 5 T, (18) 


J 
ye aye; +6,-—ri+M(1—2!)>0, t=1,2,---,T, i=1,2,...,1, (19 
TEX. Zeo, (20) 


由 于 实际 中 了 J, 了 等 的 值 均 比较 大 ， 上 述 问题 往往 是 大 规模 的 线性 整数 规划 ， 直 接 用 经 典 的 
分 枝 定 界 法 或 者 市 平面 法 很 难 在 短 时 间 内 找到 满意 的 解 。 为 有 效 求解 问题 (17)-(20)， 我 们 将 借 
鉴 处 理 大 规模 线性 规划 的 分 解 算 法 的 思想 ， 通 过 对 偶 的 方法 将 问题 的 求解 分 解 为 一 系列 小 规模 
问题 的 处 理 。 

观察 整数 线性 规划 (17)-(20)， 除 了 概率 约束 外 ， 其 余 约 束 都 因为 关联 所 有 决策 变量 而 不 具 
可 分 性 。 如 果 强 行将 其 关于 各 个 样本 分 成 了 个 子规 划 问 题 ， 则 可 能 会 导致 在 不 同情 形 下 的 最 
优 决 策 不 一 致 ， 为 克服 这 一 困难 ， 我 们 引入 类 似 随 机 规划 中 非 预 期 性 约束 (nonanticipativity 
constraints) 的 一 组 约束 为 za = 2, ii: = 1,… J, KE ot RRR HAL FRE 
量 的 值 ， 其 余 类 同 。 此 约束 要 求 任意 两 个 不 同情 景 下 的 决策 相同 。 为 减少 此 类 约束 的 个 数 ， 
我 们 可 要 求 对 每 一 个 情景 ， 都 有 一 个 约束 来 保证 其 下 的 决策 与 男 一 任意 情景 下 的 决策 相同 ， 
这 样 仅 需 要 了 个 约束 。 值 得 注意 的 是 ， 这 里 引入 约束 的 顺序 是 多 样 的 ， 如 令 任 意 情景 下 的 决 
策 均 等 于 第 一 个 情景 下 的 决策 zi = zl, i = 2,… ,TJ， 抑 或 令 任意 两 个 相 邻 情景 下 的 解 均 相 
等 zl = g”, rê = gl, i= 2, 以 下 我 们 采用 第 二 种 ， 即 每 两 个 相 邻 情景 间 进 行 约束 的 

对 应 于 上 述 引 入 的 情景 决策 ， 原 本 的 目标 函数 f(z) 就 可 以 表示 为 fat), f(z?),… ,f(z7) 的 
某 种 线性 组 合 ， 本 文 根 据 各 个 情景 发 生 的 概率 将 它们 进行 加 权 求 和 ， 即 令 


了 
Ate) =Y zf) 


有 了 上 述 的 非 预期 性 约束 ， 我 们 就 可 利用 Lagrange 松 弛 法 ， 将 所 有 非 预 期 约束 和 概率 约束 
通过 引入 Lagrange 乘 子 放 缩 至 目标 函数 中 ， 从 而 保证 剩余 约束 均 相 对 于 不 同 的 情景 ?可 分 。 该 
方法 的 优点 在 于 可 通过 一 系列 较 容 易 求解 的 小 规模 问题 的 处 理 来 逼近 原本 大 规模 、 很 难 求解 的 
问题 ， 可 大 幅度 提高 求解 效率 ， 在 较 短 时 间 内 得 到 满足 精度 要 求 的 解 。 

Lagrange 松 弛 因子 的 引入 将 原 规划 问题 转化 为 一 个 两 层 的 优化 问题 。 具 体 地 ， 





入 = [Al Aa po" Az], 和 = Dia Ata» es At S, ARs]: A= AY, A27]®, 
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内 层 需 求解 的 最 优化 问题 可 表述 为 
T 了 . | 
Vea = min fle) + SoM (Sows — ac) +> (a9 (2} - 21) + Ns 2 )) e 
t=1 i=l 


j=l i=2 
J 
st. 》 ayri tõ- rit M(L— zi) 20, t=1,2,--,7, i=1,2,. ,1 (22) 


j=l 


ze {0,1} T, z EX, i=1,.…,1. (23) 


同时 ， 外 层 的 优化 问题 需 对 Lagrange 乘 子 寻 优 ， 即 求解 下 列 优化 问题 


24 
may VLR(A) (24) 


对 于 整数 线性 规划 问题 (21)-(23)， 因 为 目标 函数 和 约束 均 可 对 i 分 解 ， 故 而 其 可 等 价 于 了 个 
子规 划 。 对 于 每 个 子规 划 ， 我 们 均 可 运用 分 枝 定 界 法 与 割 平 面 法 相 结合 的 方法 来 快速 求解 ， 如 
可 直接 利用 Xpress 等 优化 软件 。 以 下 不 再 歼 述 。 

对 于 要 求 部 分 变量 非 负 的 连续 优化 问题 (24)， 我 们 可 采用 基于 线性 搜索 的 优化 算法 。 那 么 
最 核心 的 内 容 即 在 于 如 何 确定 搜索 方向 和 选取 合适 的 搜索 步 长 。 在 本 文 的 模型 中 ， 因 原 问题 为 
离散 问题 ， 其 本 身 不 具有 梯度 ， 而 精确 的 次 梯度 又 很 难得 到 ， 故 我 们 使 用 以 下 方向 


I I I 
i i ii i i 1 Ià I 1 I 
d= [> má a,) mz- 2 TW tm — QAT, £] — L1, T3 ~ Lz, ,TI — Lz, 
i=1 ' i=l i=1 
2 1 2 1 2 1 I I-1 „I 7 一 1 了 al" 
TI ` Ti TIT To EJ EJ Ti Zi Ha — LQ 4 Hy TJ ie (25) 


这 种 方向 对 应 于 文献 中 所 谓 的 替代 次 梯度 。 对 于 一 般 线 性 整数 规划 问题 ， 文 献 22] 肯定 了 替代 
次 梯度 方向 在 Lagrange 松 弛 类 方法 中 的 有 效 性 ， 认 为 它 能 减少 对 偶 问 题 求 解 过 程 的 计算 量 ， 
这 正 是 我 们 按 上 述 方法 确定 搜索 方向 @ 的 原因 。 

虽然 我 们 选取 的 方向 可 看 作 是 对 次 梯度 方向 的 近似 ， 但 因原 问题 的 离散 性 ， 经 典 次 梯度 算 
法 中 的 步 长 规则 并 不 适用 我 们 的 问题 。 为 有 效 求解 问题 (24)， 我 们 采用 在 非 光滑 优化 中 较 常 用 
的 信赖 域 技术 B2s3,2 来 确定 搜索 步 长 。 信 和 赖 域 方法 的 思想 是 依照 每 次 迭代 对 目标 函数 值 改 进 的 
效果 来 确定 下 一 步 的 步 长 ， 若 上 一 步 欠 代 改进 明显 ， 说 明 此 方向 很 好 ， 则 加 大 步 长 以 求 更 快 下 
降 ， 若 上 一 步 迭 代 改 进 不 明显 ， 则 缩小 步 长， 以 免 在 不 好 的 方向 上 前 进 过 多 。 应 用 到 我 们 的 问 
题 ， 我 们 将 通过 实时 观察 此 次 迭代 和 上 次 迭代 之 间 改 进 的 程度 ， 这 里 选用 目标 函数 值 总 的 下 降 
量 作 为 改进 程度 的 度量 。 当 本 次 迭代 目标 值 减 去 上 次 迭代 目标 值 小 于 如 时 ， 则 增 大 步 长 ， 当 
相应 的 差 值 大 于 gr 时 ， 则 减 小 步 长 ， 当 介 于 前 述 两 种 情形 之 间 时 ， 步 长 保持 不 变 。 

至 此 ， 我 们 可 构建 如 下 确定 最 优 Lagrange 乘 子 的 迭代 算法 如 下 : 

步骤 1 YIRA Lagrange R TAO = 1 = DT VO = 0 = 4,1, j= 
1……,J， 记 相应 的 乘 子 向 量 为 X0)， 并 选取 参数 br > 1 0 < gr <1, n Mny 的 值 。 给 
定 迭 代步 数 上 限 K, 令 k=1, 5, =1- 

步骤 2 求解 对 应 于 当前 乘 子 和 *) 的 最 优化 问题 (21)-(23)， 得 最 优 解 cO KERRE Vr A). 

步骤 3 #VLr(A®P)-Vrr(A AFD) <br E) = nrs; E Ver A®)-Vrr( AFD) > 
bu, Es = ny sk- AM, Ss = sk), 

步骤 4 int (25) AR PRR dE, HEACHD = AE) + 8). dq(*)。 若 对 某 


At, AD <0, WAY = 0。 
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PRS 若 对 任意 的 1, Mj MA = ze 的 ， 则 终止 ， 当 前 的 解 向 量 O 即 为 问 
题 (17)-(20) 的 最 优 解 ; Ek > K, Aik, Ma 为 问题 (17)-(20) 的 近似 解 ， AM, Sk = 
k+1, BADR. 

对 于 联合 概率 约束 模型 (10)-(12)， 亦 可 通过 与 对 问题 (7)-(9) 类 似 的 方法 将 其 转化 为 下 列 整 
数 线性 规划 问题 


min f(z) 


I 
Tos 
s.t. dre 2o, 


azi +h&—rit+M~2')>0, t=1,2,---,T, i=1,2, I, 
j=1 


rE€EX, ze{0,1¥F. 


对 于 该 问题 ， 可 完全 同上 采用 Lagrange 松 弛 法 进一步 分 解 。 因 为 求解 算法 的 导出 完全 类 似 ， 
REMTE T. 


4 数值 例子 


为 说 明 本 文 所 建 模 型 的 合理 性 、 所 导出 样本 平均 近似 问题 的 可 解 性 与 所 设计 松弛 分 解 算 法 
的 有 效 性 与 收敛 性 ， 我 们 在 本 节 进 行 两 组 数值 实验 ， 均 以 某 个 实际 话 务 中 心 的 员工 安排 为 背 
景 ， 和 包括 我 们 在 模型 (1)-(3) 后 简单 介绍 的 公平 性 、 特 殊 班次 等 约束 。 

首先 考虑 第 一 组 实验 ， 我 们 利用 历史 话 务 量 数据 ， 预 测 出 未 来 一 月 的 基本 人 力 需 求 量 。 考 
不 到 未 来 需求 的 随机 波动 性 与 预测 的 不 精确 性 ， 我 们 假定 预测 值 存在 一 服从 期 望 值 为 0 的 正 态 
分 布 的 随机 扰动 。 也 就 是 说 ， 人 力 需 求 量 可 被 看 做 服从 均值 为 预测 人 力 需 求 量 的 正 态 分 布 。 基 
于 正 态 分 布 的 累加 分 布 函数 ， 我 们 很 容易 将 模型 (1)-(3) 转化 为 整数 线性 规划 问题 ， 并 找到 其 精 
确 最 优 解 。 

代替 上 述 通常 的 预测 方法 ， 我 们 采用 样本 平均 法 来 构建 形 如 (17)-(20) 的 规划 问题 。 有 具体 
地 ， 使 用 不 同 数量 的 样本 来 近似 人 力 需 求 量 的 分 布 ， 然 后 在 不 同样 本 集 下 求解 样本 平均 近似 问 
题 (17)-(20)， 并 将 所 得 解 与 原 问 题 的 精确 最 优 解 进 行 比较 ， 检 验 近 似 问题 的 最 优 解 是 否 随 着 样 
本 数量 的 增 大 收敛 于 原 问题 的 精确 最 优 解 。 

为 克服 简单 随机 抽样 工作 量 大 、 方 差 大 的 不 足 ， 我 们 使 用 等 距 抽 样 的 方法 来 生成 样本 
数据 集 。 这 种 方法 的 好 处 在 于 抽样 过 程 简单 。 相 较 于 抽取 了 个 样本 的 简单 随机 抽样 ， 等 距 
抽样 所 需 随 机 抽取 的 次 数 仅 为 1， 这 节约 了 很 多 工作 量 。 我 们 分 别 选取 20 个 样本 与 50 个 样 
本 进行 实验 ， 同 时 与 原 问题 进行 对 比 。 三 组 实验 采用 相同 的 置信 度 与 偏差 容忍 度 ， 并 对 任 
Mt=1,---,T, Ma, = 0.98, 6 = 一 2 + 0.05 * r JETER, RIER 1 中 列 出 实际 中 各 
天 所 需 的 最 小 人 力 。 


表 1: 样本 近似 问题 最 优 值 与 原始 问题 最 优 值 的 比较 








天 1 2 3 4 5 6 7 8 9-31 ”合计 
原始 问题 84 84 77 70 68 66 68 68 wee 2209 
20 个 样本 平均 近似 87 88 81 73 7 69 71 7 T 2309 


50 个 样本 平均 近似 85 85 78 71 68 66 68 68 eri 2222 
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为 对 比 说 明 不 同 结果 ， 我 们 将 不 同 天 原 问 题 精确 的 最 优 值 、20 和 50 个 样本 下 样本 平均 问题 
的 最 优 值 ， 即 每 天 所 需 最 小 人 力 绘制 于 图 1 中 。 

从 图 1 中 不 难看 出 ， 随 着 样本 量 的 增 大 ， 样 本 平均 近似 问题 的 最 优 值 逐 渐 逼 近 原 始 问题 的 
最 优 值 ， 尤 其 是 在 50 个 样本 时 ， 差 别 就 已 经 很 小 了 。 








l 
ST 


图 1: 样本 平均 近似 问题 最 优 值 与 原 问题 最 优 值 的 比较 


作为 用 户 可 调节 的 选项 ， 偏 差 容忍 度 决 定 了 概率 约束 在 限制 可 行 解 上 所 起 作用 的 大 小 ， 从 
而 影响 到 问题 最 优 值 的 大 小 。 为 此 ， 我 们 考虑 另 一 组 实验 ， 研 究 不 同 偏差 容忍 程度 对 问题 最 
优 解 与 最 优 值 的 影响 。 我 们 考虑 两 种 形式 的 偏差 : 绝对 偏差 和 相对 偏差 。 不 难 理解 ， 绝 对 ( 相 
对 ) 偏差 越 少 ， 对 应 问题 的 最 优 值 越 小 ， 而 相对 偏差 方式 较 绝对 偏差 方式 ， 对 应 的 最 优 值 一 般 
较 大 。 

为 验证 上 述 结论 ， 我 们 设计 了 四 组 实验 ， 它 们 分 别 是 ， 较 大 的 绝对 偏差 biast = 1, t= 
Lie, Ts 较 小 的 绝对 偏差 biast = 0.5, t= 1,… Ts 较 大 的 相对 偏差 biast = 0.05 «74, t= 
1,.… ,了 和 较 小 的 相对 偏差 Dias, = 0.08 * rt，t = 1,… ,人 。 我 们 再 统一 令 6, = -2 + bias, t = 
1 ,了 7 了 ， a, 二 0.98, t= 1,… T, EREI = 20。 为 考虑 概率 约束 不 起 作用 时 问题 最 优 值 的 
变化 ， 我 们 特 增加 一 次 对 照 实验 ， 它 对 应 于 bias, = -~co， 此 时 所 得 最 优 值 应 为 上 述 情形 下 最 
优 值 的 下 界 。 上 述 各 种 情况 下 所 得 最 优 值 见 表 2。 


表 2: 不 同 偏差 容忍 程度 下 问题 最 优 值 的 比较 











天 1 2 3 4 5 6 7 8 931 合计 

bias, = 一 oo (无 约束 问题 ) 78 78 72 65 63 61 63 63 ... 2055 
bias, = 0.1 8&2 81 75 68 66 64 66 66 =- 2215 

bias, = 0.5 80 81 75 68 66 64 66 66 ... 2140 

bias; = 0.05 * rt 87 88 81 73 71 69 71 71 ... 2309 
bias; = 0.03 * re 82 84 78 70 68 67 69 69 > 2218 





为 进一步 了 解 整个 月 内 最 优 值 的 变化 ， 我 们 将 上 述 五 组 解 的 结果 在 图 2 中 予以 展示 。 

从 表 2 与 图 2 中 可 以 看 出 ， 无 论 是 绝对 偏差 还 是 相对 偏差 ， 偏 差 值 较 大 时 问题 的 最 优 值 均 较 
大 ， 较 小 时 间 题 的 最 优 值 较 小 ， 这 与 前 述 的 理论 结论 是 一 致 的 。 而 绝对 偏差 情形 与 相对 偏差 情 
形 相 比较 ， 绝 对 偏差 下 对 应 的 最 优 值 要 小 于 相对 偏差 下 对 应 的 最 优 值 。 
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| 一 -一 biasi=-inf 








ae ks bias=1 


ee bias =0.5 
~— bias =0.05*r, 
bias =0.03"r, 











B2: 不 同 偏差 容忍 程 度 下 问题 最 优 值 的 比较 


结论 


相 较 已 有 文献 ， 本 文 所 建立 的 带 有 概率 约束 条 件 的 智能 排 班 模型 更 精确 地 描述 了 人 力 需求 


限制 的 要 求 ， 从 大 概率 角度 有 效 地 限制 了 人 力 不 足 ， 更 符合 客户 期 望 。 而 我 们 基于 样本 平均 近 
似 方法 、Lagrange 松弛 法 、 分 解 算 法 等 思想 所 设计 的 堪 套 分 解 算 法 也 较 已 有 算法 更 为 有 效 。 然 
而 ， 在 样本 平均 近似 方法 中 ， 样 本 的 选取 对 结果 有 着 较 大 的 影响 ， 故 而 如 何 寻找 一 种 更 稳定 、 
有 效 的 抽样 方法 仍 待 更 深入 的 研究 。 
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Chance Constrained Programming Models for Intelligent Scheduling 
Problems and Their Efficient Solution 
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(1- Department of Scientific Computing and Applied Software, Faculty of Science, 
Xi’an Jiaotong University, Xi’an 710049; 2- Jinhe Center for Economic Research, 
Xi’an Jiaotong University, Xi’an 710049) 


Abstract: Due to their extensive application areas and solution difficulty, different forms of intelli- 
gent scheduling problems have now become a research hot in operations research. The current models 
can not properly reflect practical characteristics of labor requirement constraints and the relevant al- 
gorithms lack necessary efficiency. To overcome these shortcomings, two forms of chance constrained 
programming models are proposed for comprehensive scheduling problems to flexibly describe the labor 
requirement constraints. By using the sample average approximation technique and the Big-M method, 
we transform the chance constrained programs into large-scale integer programming problems. An ef- 
ficient nested decomposition algorithm is designed by simultaneously utilizing the Lagrange relaxation 
technique, the decomposition method and the trust region scheme. Finally, some numerical results 
illustrate the reasonability of the proposed models and the efficiency of the proposed algorithm. 
Keywords: intelligent scheduling; labor requirements; chance-constrained; sample average approxi- 
mation; the big-M method; Lagrangian relaxation; trust region 





Received: 12 Aug 2010. Accepted: 18 Oct 2010. 
Foundation item: The National Natural Science Foundation of China (70971109). 


